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1. Introduction 

Central ultrarelativistic collisions at RHIC and LHC with more then 3000 
produced particles present remarkable opportunity to analyse event-by-event fluctu- 
ations of hadronic observables. Such single events analysis with large statistics can 
reveal new physical phenomena usually hidden when averages over a large statistical 
sample of events are made [1,2]. The numbers of particles produced in relativistic 
heavy-ion collisions can differ dramatically from collision to collision due to the 
variation of impact parameter (centrality dependence), energy deposition (leading 
particle effect), and other dynamical effects [2,3]. The fluctuations can also be in- 
fluenced by novel phenomena such as disoriented chiral condensate (DCC) [4,5]. 

The best probes of such novel dynamics are fluctuations of conserved quanti- 
ties, because conservation laws limit the degree to which final-state scattering can 
dissipate. Even globally conserved quantities such as energy, net charge, isospin, 
baryon number and strangeness can fluctuate when measured, e.g., in a limited 
phase space region. 

Several methods have been proposed [3,7] to distinguish between statistical 
and dynamic fluctuations . In high energy hadronic and heavy-ion collisions the 
correlations of particle production are usually analysed and described in terms of 
"short range correlations" (SRC) and "long range correlations" (LRC). The SRC 
are mainly due to clustering of outgoing particles. These clusters may be hadronic 
resonances which decay into a few particles or minijets [5,6,8]. The LRC dominate 
in high multiplicity events and are mainly due to global conservation laws. 

The old puzzle in cosmic-ray observations is the existence of few exotic events 
characterized by an anomalously large number of charged pions in comparison with 
the number of neutral pions, the Centauros [9], indicating that there should exist a 
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strong negative long-range correlation between two types of the pions. 

Such long-range correlations are possible if pions are produced semiclassicaly 
and constrained by global conservation of isospin [10-15]. 

Although the actual dynamical mechanism of the production of a classical 
pion field in the course of a high-energy collision is not known, there exist numer- 
ous interesting recent theoretical attemps to explain Centauros either as different 
types of isospin fluctuations due to the formation of a disoriented chiral condensate 
(DCC) [16-20], or as multiparticle Bose- Einstein correlations (BEC) [21 ], or as 
the formation of a strange quark matter (SQM) [22]. Among the most interesting 
speculations is the idea of DCC that localized regions of misaligned chiral vacuum 
might occur during the ultrahigh-energy hadronic and heavy-ion collisions when the 
chiral symmetry is restored at high temperatures. These regions, if produced, would 
behave as a pion laser, relaxing to the ground state by coherent pion emission [ ]. 
It is generally accepted that the fluctuation of the ratio of neutral to charged pions 
of the Centauro type could be a sign of the DCC formation provided that a single 
large domain is formed, containing a large number of low Pt pions. Since the pions 
formed in the DCC are essentially classical they form a quantum superposition of 
coherent states with different orientation in isospin space. If all the pions in the 
domain are pointing in the same isospin direction and the condensate state is a pure 
isoscalar then the formation of DCC leads to large event-by-event fluctuations in 
the ratio f — no/n oi the number of ttq's in the DCC divided by the total number 
of pions produced in an event. The probability distribution of / inside the DCC 
domain is [13,16,17,19] 

Poccif) = ^ (1) 
There are a variety of proposed mechanisms other than DCC which also lead to the 
distribution (1) [23,24,25,26]. The distribution Poccif) different from the generic 
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binomial-distribution expected in normal events which assumes equal probability for 
production of 7r+, 7r_ and ttq pions. The emission of charged and neutral pions is 
then uncorrelated and the distribution 



in the limit as n — > oo, no — >^ oo with / fixed, approaches a delta function at 



The possibility of observing the DCC type fluctuations critically depends on 
the size and the energy content of the DCC domain. If the domain is of the pion size, 
the effect of DCC is too small to be observed experimentaly. The early accelerator 
searches for Centauros and DCC at CERN [27-30] and at Fermilab [23,31] were thus 
unsuccesful. With the RHIC facility at BNL now, there is a possibihty to consider 
event-by-evevnt fluctuations of other speciflc hadronic observables that might be 
more informative about the signal of DCC. There are various factors that may affect 
the observation of the DCC signal [32]. Experimental signals such as the isospin 
fluctuations, the strong relative enhancement of the number of low pr pions, and 
the suppresion of HBS correlations may provide a robust signal of DCC [33]. 

The space-time scenario of the formation and decay of the DCC is usually 
studied within one of the simplified versions of the chiral effective Lagrangians, either 
the linear or nonlinear sigma model [20]. However, it should be emphasized that 
the use of cr- models, be they linear or nonlinear, is only a rough approximation to 
the true dynamics, because the couplings of pion and sigma fields to the constituent 
quarks may be large and their effect should not be ignored. 

In this paper, we present results of an event-by-event analysis of charge-neutral 
pion fluctuations as a function of the p/tt ratio in pp-coUisions. Following the ap- 
proach of our earlier papers [24,34], we consider in Section 2. the leading-particle 





/ = 1/3. 
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effect as a source of a classical pion field in the impact parameter space which is 
related to the semiclassical solutions of the equation of motion of the nonlinear a- 
model coupled to quark degrees of freedom. In order to faciliate the analysis of 
charged-neutral pion correlations, we also derive the corresponding pion-generating 
function. The quantum equations of the nonlinear cr-model coupled to quarks are 
discussed in Section 3., and shown to lead to a coherent state description of the 
pions emerging from the DCC. A multipion exchange potential between two quarks 
is derived as well as a corresponding correlator for multipion production. Results of 
our investigation are summarized in the Section 4. Our general conclusion is that 
within the nonlinear a model the large isospin fluctuations depend strongly on the 
value of the p/tt ratio which fluctuate from event to event. 

2. Pion production from a classical source 

At high energies most of the pions are produced in the nearly baryon-free 
central region. The energy available for the hadron production is 

Ehad — -^Vs — Eieading (3) 

which at fixed total cm. energy \/s varies from event to event. The n-pion con- 
tribution to the s-channel unitarity can be written as an integral over the relative 
impact parameter h of the two incident leading particles: 

1 /■ " 

F„(s) ^- d%X{ dq, I 7;(s, b;l...n)\', (4) 
where dq = d'^qTdy /2{2t:)^. The normalization is such that 

F„(s) = sa„(s), 

oo 

CTinelis) = ^(T„(s). (5) 

n=l 

If the isospin of the incoming leading particles is II3 , then the initial-state vector 

A — * 

of the pion field is S{s, b) \ II3), where | 7/3) denotes a vacuum state with no pions. 
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The n-pion production amplitude is 

iT^is, b;qi... Qn) = 2s(/7^; gi . . . g„ | S{s, b) \ Ih), (6) 

where I'l'^ denotes isostate of the outgoing leading particles. 

The basic assumption of the independent emission of pion-clusters, in b-space 
is the factorization property of the scattering amplitude. The factorization of T„ 
follows if the quantum field of the pion-cluster satisfies the equation of motion of 
the form 

{a + ml)7rc{s,b;x) = fc{s,b;x), (7) 

— * 

for an isovector cluster and similar one for an isoscalar cluster. Here, jc is a classical 
source of the cluster which decay into c = 1,2,... pions outside the region of strong 
interactions (the final-state interaction between pions being neglected). Clusters 
decaying into two or more pions simulate a short-range correlation between pions. 
They need not be well-defined resonances. 

If the conservation of isospin is a global property of the colliding system, then 
jc{s, b; x) is of the form 

jc{s,b;x) ^ jc{s,b;x)n, (8) 

where n is a fixed unit vector in isospace independent of x, s and b . The global con- 
servation of isospin thus introduces the long-range correlation between the emitted 
pions. 

The S matrix following from such a classical source is still an operator in the 

A — * 

space of pions. Inclusion of isospin requires S{s, b) to be also a matrix in the isospace 
of the leading particles. 

The coherent production of pion-clusters is described by the following S matrix: 
S{s, b)^ J (ffi I n)D{s, b){n\, (9) 
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where | n) represents the isospin-state vector of the two- leading-particle system. 
The quantity D{s, h) is the unitary coherent-state displacement operator defined in 
our case as 

D{s, h) = exp[o){s, h) - a(s, h)] (10) 

where 

h)^Y.I HUs. b; q)fiaj{q) + Z{s, b; q)al{q)], (11) 

and aj{q) and aj(g) are the creation operators of a cluster of type c, respectively. 

The isospin {I'jI's) of the outgoing leading particle system varies from event 
to event.lt is produced with the probability ^vj'^, and we can sum over all (Z'/s) to 
obtain the probability distribution of producing n+Tr"*", n_7r-, and no7r° from a given 
isospin state: 

Pii3{n+n_no)Nii^ = / d%dqidq2 . . . dqn \ {I'l'^n+njiQ \ S{s,b) \ II3) 

^ ^3 

n = + n_ + no (12) 
where Nu^ is the corresponding normalization factor determined by 

E{n}^//3(^+>^->^0) = 1- 

This is our basic relation for calculating various pion-multiplicity distributions, 
pion multiplicities, and pion correlations between definite charge combinations. In 
general, the probability Pii^{n+n_no) would depend on {I'l'^) dynamically. The 
final-leading-particles tend to favor the {I' ,1'^) ~ {I,h) case. However, if the leading 
particles are colliding nuclei, an almost equal probability for various (/', Jg) seems 
reasonable aproximation owing to the large number of possible leading isobars in 
the final state. Then we can sum over all {I'l'^) using the group theory alone. 

Recent studies of heavy-ion collisions at the partonic level [26] argue that the 
central region is mainly dominated by gluon jets. The valence quarks of the incom- 



7 



ing particles which escape from the interaction region form the outgoing leading 
particle system. Since gluon's isospin is zero, it is very likely that total isospin of 
the produced pions in the central region is also zero. This picture is certainly true 
if the central region is free from valence quarks, the situation expected to appear at 
the extremelly high collision energies. 

Let us now assume that pions are produced both singly and through isovector 
clusters of the p type [24]. In this case, the most appropriate tool for studying 
various pion correlations is the generating function Gii^{z,n_): 

Gn,{z,n_)= ^ P,,3(n+,n_,no)/^°, (13) 

no,n+ 

from which we can calculate, for example 

{no)n_ = ^lnGii^{l,n_), (14) 



and 



fln_ = ^lnGns{l,n_), (15) 



1 

Pn.{no) = -^^X^uMn_). (16) 



n 



The form of this generating function is 

G„,(.,n.) = (/+ /' I F/'W 1^ i(£l£)!!e-^(-.-), (17) 

2 (7 + 73)! J -I UqI 

where 

2A{z, x) = (1 - x^)n^ + z{l - x^)np + 2x% (18) 

and 

2B{z, x) = n^{l + - 2zx^) + np{2 - z{l - x^)). (19) 

Here n.^ denotes the average number of directly produced pions, and rip denotes 
the average number of p-type clusters which decay into two short-range correlated 
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pions. The function Pj^{x) denotes the associate Legendre pohnomial. Note that 



A{l,x) = B{l,x). 



The total number of emitted pions is 



n — + 2np. 



(20) 



The behavior of P(no) = -Poo(^o) for n = 50 and different combinations of 
(rr,r,^p) is shown in Fig. l.The behaviour of the no-dispersion 



for a given number of negative pions and different pairs of (TTtt, Tip) in the case of 
I — — 1 is shown in Fig. 2.. Note that /2,n_ is a sensitive quatity of the pairing 
properties of the pions. 

We see that DCC-type behavior is expected only for 7^ and rip — 0. that is 
in events without p-resonances. Recent estimate of the ratio of p- mesons to pions, 
at accelerator energies, is rip = O.lOn^r [6]- 

3. DCC from a quantum nonlinear sigma model 

In this section we establish the relationship between the quantum nonlinear 
sigma model coupled to quarks and our cohetrent state eikonal model description of 
the pions emerging from the DCC. 

As is well known, the Lagrangian for QCD with two light up and down quarks 
has an approximate global SU{2)l x SU{2)r symmetry, which at low temperatures, 
is spontaneously broken to SU{2)v by a nonzero value of the quark condensate 
(QlQr), which is regarded as an order parameter of the system. This order parameter 
can be represented as a four-component vector 4> = (c, ^) buildt from the quark 
densities. The chiral symmetry then corresponds to 0(4) rotations in internal space. 



D{no)l_ = h,n_ - {no) 



n. 



(21) 
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The true vacuum of the theory is defined as (0) = {{a), 6), with (a) ^ 0. 
In QCD the spontaneous symmetry breakdown leads to nearly massless Goldstone 
bosons (the pions) and gives the constituent-quark mass. At low energies and large 
distances ( momentum scale smaller than \GeV ) the dynamics of QCD is described 
by an effective Lagrangian containing the cr, tt fields and constituent quarks. 

In the DCC dynamics we distinguish three stages: formation, evolution and 
decay stage. In the conventional approach [19] one starts with a chirally symmetric 
phase at TT^ and DCC formation happens as T, due to a rapid expansion or cooling, 
drops below T^. spontaneously breaking the chiral symmetry. 

The evolutionary stage of the DCC is usually described by the classical chi- 
ral dynamics based on the cr-model, mostly the linear a- model. We consider the 
nonlinear cr-model coupled to quarks at zero temperature [35] which is expected to 
describe the late stage of the DCC evolution. 

The Lagrangian for the nonlinear cr-model coupled to quarks is 

L = ^Trid^U^d'^U) + q{nd)q - gUqUq (22) 

where 

U = exp(i75^) (23) 

We parametrize the pion field in the following form 

7f(x) = Un{x)e{x) (24) 

where n{x) is an unit vector which determines the isospin orientation of the pion 
field, obeying — 1. 

The Euler-Lagrange equations of motion for 9 and n are 

aO - sin 6 cos Od^n ■ d^n = -i-^n ■ {Qf'^^Q) 
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df.isiii'en X d^n) = -i^n x (Qr-f^Q) sin 9 (25) 

(26) 

where Q denotes the constituent quark defined by 

Q = e'^'Mq (27) 

and tuq = qJt^ is the constituent quark mass. 
We treat 

-i^gT75g-J(x) (28) 

as a given classical external source of pions and identify it with the source function 
ji (s, 6; x) in our eikonal model. For the class of solutions that can be rotated into 
a uniform one, n{x) — n, known as the Anselm-class of solutions [16], the solutions 
with constant n can be realized if the source points to a certain fixed direction n in 
the isospace: 

j{x)^j{x)n (29) 
Then the equation of motion for the pion field reduces to 

Deix)^j{x) (30) 

with 

7r{x) = UO{x)n (31) 

This relationship offers us the possibility to study the importance of various 
quark sources in the DCC formation. 

In quantum chiral field theory of the nonlinear cr-model the role of a strong 
coupling of gTT-interaction has the quantity mq / Jt^. Since the Lagrangian of the non- 
linear (T-model is nonpolynomial, a suitable renormalization procedure and operator 
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normal ordering should be formulated [29] . Let 

$ =: - 1 : (32) 

denote the chiral super field of the pion. The Lagrangian describing the interaction 
of this super chiral field with quarks is now 

L = ^Tridf.^d"^) + q{i^d - mQ)q - ^q^q. (33) 

The chiral super propagator of the field $ is defined by 

A^{x) = (r($(x)$(0))) (34) 

Its explicite form, after number of algebraic manipulations, is 

A^{x) = {1 ® 1 + ^(75t) (75t)9a}9a[Ac/i(^)] + 1 (35) 

Jit 

where 

(T(7r,(x)7r,(0))) = 5,,A(a;). (36) 

and 



A{x) = 



47r^ — ie 

The multipion exchange potential between two quarks is related to the Fourier trans- 
form of A$(x) in the following way: 

u{p[)u{p2)A^{q)u{pi)u{p2) = u;[uj2V^{q)uJiUJ2 (37) 

where 



uip) 



lep + m uj^ 
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and 



ap 
ep + m 



-UJ 



with Oi)*uj — 1. 

For studying the multipion production on quarks, it is necessary to find the 
correlator 

Anm(9l, ?2, Qn] h: k: ...Jni\x) = {n\ T($(a;)$(0)) | TJl) (38) 

This is an involved task and has not been done yet for arbitrary n and m. 

In the case of soft chiral pion bremsstralung [37-40], in which every incoming 
and outgoing quark line is replaced by 

— * — * 

TTT 

q — ^ exp{i-f5—)q (39) 

the distribution of neutral pions for large n is of the form 1 / ^n^n which is typical 
for coherent pion production without invoking the notion of DCC formation. 

Therefore, we conclude that nP(no | n) ~ ^Jn/riQ behaviour is not always a 
definite signature of DCC formation. 

4- Conclusion 

The results of the present analysis have shown the experimental observation 
of DCC is strongly affected by the p/tt ratio. In particular, we have found that: 

• Within the framework of an unitary eikonal model with factorization, energy 
conservation and global conservation of isospin the DCC-type fluctuation of 
the neutral pion fraction (/) could be observed if the p/tt ratio is small in an 
event; 
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• The DCC effect also depends on isospin of the initial-leading- particle system; 

• The coherent production of p-type clusters supresses the DCC-type behaviour; 

• The factorization property of the scattering amplitude in the impact parameter 
space may be related to the isospin-uniform solutions of the quantum nonlinear 
cr-model coupled to quarks; 

• The multipion exchange potential between two quarks can be found; 

• The soft chiral pion bremsstralung also leads to anomalously large ratio of 
neutral to charged pions. 

Acknowledgment This work was supported by the Ministry of Science 
of Croatia under Contract No. 0098004. 
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Figure captions : 

Fig.l. The curves represent P{no) for different combinations of (n7r,np), the 
average number of singly produced pions and the average number of p-type clusters, 
respectively. 

Fig. 2. The dispersion for two neutral pions as a function of the number 
of negative pions for 7 = 73 = 1. The curves represent different combinations of 
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